Nuclear signatures in high-harmonic generation from laser-driven muonic atoms 
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High-harmonic generation from muonic atoms exposed to intense laser fields is considered. Our 
particular interest lies in effects arising from the finite nuclear mass and size. We numerically perform 
a fully quantum mechanical treatment of the muon-nucleus dynamics by employing modified soft- 
core and hard-core potentials. It is shown that the position of the high-energy cutoff of the harmonic 
spectrum depends on the nuclear mass, while the height of the spectral plateau is sensitive to the 
nuclear radius. We also demonstrate that 7-ray harmonics can be generated from muonic atoms in 
ultrastrong VUV fields, which have potential to induce photo-nuclear reactions. 
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Muonic atoms have ever since been useful tools for nu- 
clear spectroscopy employing atomic-physics techniques. 
Due to the large muon mass and the correspondingly 
small Bohr radius, the muonic wave function has a large 
overlap with the nucleus and thus effectively probes its 
structural features, such as finite size, deformation, sur- 
face thickness, or polarization Muonic atoms also 
play an important role as catalysts for nuclear fusion Q • 
On the other hand, there is the steadily growing field 
of laser-matter interactions. While traditionally deal- 
ing with the atomic or molecular response to an exter- 
nal light wave, it has recently widened to include also, 
among other new directions, laser-nuclear physics 
With laser- heated clusters and plasmas, nuclear fission 
Q and fusion [B| have been realized as well as neutron 
production in (p,n) or (d,n) reactions By virtue of 
the interaction with an applied laser field, high-precision 
experiments have been able to monitor ultrafast vibra- 
tions of the nuclei in diatomic molecules 0- The rapid 
progress in laser technology even opens prospects for nu- 
clear quantum optics via direct laser-nucleus coupling Q . 

Against this background, the combination of muonic 
atoms with laser fields in order to dynamically probe nu- 
clear properties by the laser-driven muon appears very 
promising. The information on the nucleus, dynami- 
cally revealed by laser assistance, can in principle comple- 
ment the knowledge gained from field-free spectroscopy 
of muonic atoms. Modifications of muon-catalyzed fu- 
sion from a muonic molecule in the presence of a 
superintense laser field have recently been addressed Q. 
We note that on the typical time scale of strong laser 
pulses (r ~ fs— ns) muonic atoms and molecules may be 
considered as stable since the muon lifetime is 2.2 /j,s. 

In this Letter, high-harmonic generation (HHG) [13] 
from muonic hydrogen isotopes in strong laser fields is 
studied. We demonstrate nuclear signatures arising from 
the finite nuclear mass and size in the radiation emitted 
by the muonic charge cloud, which oscillates across the 
nucleus under the influence of the laser field. The rel- 
ative motion of the muon and the nucleus is taken into 
account. In the ground state of muonic hydrogen, the 



muon is bound by 2.5 keV and experiences an electric 
field strength of 1.8 x 10 14 V/cm corresponding to the 
intensity 4.2 x 10 25 W/cm 2 . It is of interest to com- 
pare these numbers with the parameters of the most 
advanced present-day and near- future laser sources. In 
the optical and near-infrared frequency range (fku ~ 1 
eV), the next generation of high-power lasers will reach 
the intensity level of 10 23 W/cm 2 0. At VUV fre- 
quencies (Hut ~ 10—100 eV) a record intensity of al- 
most 10 16 W/cm 2 has recently been achieved at the 
FLASH facility (DESY, Germany) with a free-electron 
laser 12]. Near-future upgrades of such machines are 
envisaged to generate x-ray beams (hw ~ 1 — 10 keV) 
with peak intensities close to 10 20 W/cm 2 . There are 
also efforts to produce ultrashort, high-frequency radia- 
tion (fiu) ~ 10 — 1000 eV) from plasma surface harmon- 
ics where, due to high conversion efficiency, considerably 
higher intensities might be reachable [13j. With these 
high- intensity and/or high-frequency sources of coherent 
radiation it will become possible to influence the dynam- 
ics of (light) muonic atoms. As an example, we investi- 
gate high-harmonic emission from muonic hydrogen and 
deuterium in intense VUV laser fields. Our numerical 
calculations reveal pronounced differences in the HHG 
spectrum when changing the nuclear mass and size. 

We consider muonic quantum dynamics in the non- 
relativistic regime described by the time-dependent 
Schrodinger equation (TDSE). Due to the large muon 
mass, the nucleus cannot be treated as an infinitely 
heavy particle and its motion must be taken into account. 
Within the dipole approximation for the laser field, the 
two-particle TDSE separates into center-of-mass and rel- 
ative motion. By employing a one-dimensional (ID) 
hydrogen-like model atom (to be justified below), the 
center of mass evolves in time according to 
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with the muonic, nuclear, total and reduced masses m^, 
m n , M — nip + m n and m r = m^rrin/M, respectively, 
center-of-mass and relative coordinates X and x, cor- 
responding momentum operators P = —ihd/dX and 
p = —ihd/dx, nuclear charge number Z , nuclear poten- 
tial V(x) and laser electric field E(t). Equation (Q} is the 
non-relativistic Volkov equation for a particle of charge 
(Z — l)e and mass M. In the present paper we restrict 
ourselves to the case Z = 1 as the required laser inten- 
sities are the smallest then, while the effects of interest 
are expected to be appreciable. Note that throughout 
the nuclear chart the largest relative difference among 
isotopes, both in mass and charge radius, exists between 
hydrogen nuclei. In this situation, Eq. (1) describes free 
motion of the center-of-mass coordinate which does not 
generate radiation and can be ignored. The laser field 
solely couples to the relative coordinate. For Z = 1, 
Eq. ^ is reduced to the usual Schrodinger equation for 
a single particle of charge — e and mass m r in the com- 
bined fields of a laser and a nucleus. By applying the 
standard scaling transformation, one can recast Eq. ([2]) 
into a form that describes an ordinary (i.e. electronic) 
hydrogen atom in a scaled laser field. The scaled laser 
frequency and field strength are given by Q 
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where p = m e /m ri with the electron mass m e . That is, 
a muonic hydrogen atom in a laser field with parameters 
E, uj behaves like an ordinary hydrogen atom in a field 
with E',u>' given by Eq. ([3]), provided that the Coulomb 
potential V(x) arises from a pointlike nucleus. The scal- 
ing procedure, however, does not account for nuclear pa- 
rameters like, e.g., the finite nuclear size. Evidently, 
when the transition from a muonic to an ordinary hy- 
drogen atom is performed, the nuclear radius is not to 
be scaled but remains fixed. Hence, for atomic systems 
where nuclear properties play a role, not all physical in- 
formation can be obtained via scaling. Below we show 
results which display the influence of the finite nuclear 
size on the process of HHG. 

We solve the TDSE in Eq. © numerically by making 
use of the Crank-Nicolson time-propagation scheme. The 
laser field is always chosen as a 5-cycle pulse of trape- 
zoidal envelope comprising one cycle for linear turn-on 
and turn-off each. The harmonic spectrum is obtained by 
taking the Fourier transform of the dipole acceleration. 
With regard to the ID approximation in Eqs. (|1I2|) we 
note that, in general, ID models are known to retain the 
essential physical features of nonrelativistic laser-atom 



interaction for linearly polarized fields [14(; therefore, 
these models are widely used (l5| . Furthermore, in our 
case the ID numerics still are a non-trivial task because 
of the fine grid spacing required to resolve the nuclear ex- 
tension and the non-standard laser parameters employed. 
As regards ordinary atoms, the latter would correspond 
to intense fields in the far-infrared. Eventually, we stress 
that the goal of this paper is to reveal relative differences 
between physical observables which typically are less sen- 
sitive to model assumptions than absolute numbers. The 
main shortcoming of the ID approach consists in neglect- 
ing the muon's transversal wave-packet spreading which 
can substantially reduce the total harmonic yield [lij ]. 

As the conversion efficiency of HHG is rather low, it 
generally is desirable to maximize the harmonic signal 
strength. In the present case, the optimization is of par- 
ticular importance as the target density is low. Efficient 
HHG requires both efficient ionization and recombina- 
tion. The former is guaranteed if the laser peak field 
strength lies just below the border of over-barrier ioniza- 
tion (OBI) where the Coulomb barrier is suppressed all 
the way to the bound energy level by the laser field, i.e. 
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corresponding to the intensity / OBI « 1.6 x 10 23 W/cm 2 . 
Here, a ~ 1/137 denotes the fine-structure constant. 
Efficient recombination is ensured if the magnetic drift 
along the laser propagation direction is negligibly small, 
which limits the relativistic parameter to [171 ] 
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with the atomic ionization potential I p . We note that 
condition ([5]) also implies applicability of the dipole ap- 
proximation. The lowest frequency which simultaneously 
satisfies Eqs. ^ and ([5]) lies in the VUV range and corre- 
sponds to hui ks 27 eV. At this value, the harmonic radi- 
ation has a remarkable maximum energy of E c 0.55 
MeV. Hence, muonic atoms are also promising candi- 
dates for the generation of hard x-rays which might be 
employed to trigger photo-nuclear reactions. 

To express the possible effects due to the finite nuclear 
mass appearing in HHG spectra of muonic hydro gen iso- 
topes, we introduce the usual soft-core potential 
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FIG. 1: (color online). HHG spectra calculated with the soft- 
core potential ((6]). The black and grey (red) lines represent 
the spectrum for muonic hydrogen and deuterium, respec- 
tively. The laser parameters are / = 1.05 x 10 23 W/cm 2 and 
hu) = 59 eV. 



in Fig. [I] cannot be attained by present-day technology 
in the laboratory. One might, however, employ a rela- 
tivistic beam of muonic atoms counter-propagating the 
laser beam [l9j . Then, via the relativistic Doppler shift, 
the parameters of Fig.[T]can be reached by means of an in- 
tense optical laser pulse of hu> « 1.5 eV and I s=s 6.5 x 10 19 
W/cm 2 at an atomic Lorentz factor of 7 ss 20. In their 
rest frame, the atoms would then experience the laser 
parameters of Fig. [I] For muonic hydrogen atoms car- 
rying zero charge, the acceleration cannot directly be 
accomplished, though. But one could accelerate muons 
and protons/deuterons separately and let them combine 
in a merged beams setup. A similar way of producing 
positronium atoms at 7 ~ 20 via acceleration of negative 
positronium ions has been proposed recently [20] . 
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appropriately scaled to the muonic system. Here, we con- 
sider only the variation of nuclear mass assuming the nu- 
cleus as being point-like. Figure [1] shows the HHG spec- 
tra for muonic hydrogen and deuterium in an ultra-strong 
VUV field of the frequency Hcu — 59 eV and the intensity 
I = 1.05 x 10 23 W/cm 2 . The spectrum of hydrogen is en- 
hanced and further extends by about 60 harmonics. The 
reason for the spectral enhancement is that the effective 
particle of the proton-muon system experiences stronger 
acceleration since it has a smaller reduced mass. The 
difference in the cutoff position E c = I p + 3.17U P [l(| is 
also caused by the reduced mass which enters the pon- 
deromotive energy 
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For the chosen laser parameters we obtain U p ^> I p so 
that the smaller reduced mass leads to a higher cutoff 
energy. Note that in a situation where U p <C I p , the 
relative order of the cutoff positions would be reversed. 
Furthermore, we point out that the pondermotive energy 
in Eq. ([7]) can be understood as the sum of the pondero- 
motive energies of muon and nucleus. As schematically 
depicted in Fig. [21 both particles are driven into opposite 
directions by the laser field and when they recollide their 
kinetic energies add up. In this picture the higher cut- 
off energy of the hydrogen atom arises from the larger 
ponderomotive energy of the proton as compared to the 
deuteron. Note that the phenomenon that both atomic 
constituents are driven into opposite directions is most 
pronounced for positronium, the bound state of an elec- 
tron and positron [lj|. We also observe a few kinks in 
the spectra of Fig.[TJ This multiplateau structure is typ- 
ical of the interaction regime where £ ~ 0.1 (in our case 
£ w 0.04) 0. 

At the chosen frequency, the laser intensity assumed 



Muonic Deuterium 

FIG. 2: Schematic diagram of the motion of the nucleus (on 
the left) and the muon (on the right) in muonic hydrogen and 
deuterium in the presence of a laser field polarized in horizon- 
tal direction. The oscillation amplitudes are not scaled. 

To express the possible effects due to the finite nuclear 
size appearing in HHG spectra of laser-driven muonic 
hydrogen isotopes, we need a potential which explicitly 
takes the nuclear radius into account. We make use of the 
nuclear drop model, considering the nucleus as a sphere 
of uniform charge density within the nuclear radius i?, 
and employ the potential: 
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We note that the lowest lying atomic state of this poten- 
tial is unphysical as its binding energy becomes infinite 
in the limit R — > |2l |. Therefore, in accordance with 
common practice [2lL I22L l23l |. we start our calculation 
from the first excited state which has the correct hydro- 
genic binding energy. Figure [3] shows the HHG spectra 
obtained from the potential |(5J) for muonic hydrogen and 
deuterium. The proton and deuteron charge radii are 
R p m 0.875 fm and R d « 2.139 fm 0]. In order to 
compensate for the nuclear mass effect, we apply accord- 
ingly scaled frequencies and intensities with uj oc m r and 
I oc m\ [cf. Eq. J5])] so that the cutoff positions of both 
spectra coincide. We note that the overall shape of the 
harmonic response differs from that shown in Fig. [I] due 
to the different potentials employed. In particular, a dip 
can be seen at low harmonics (nfko ~ I p ) followed by a 
rising plateau region. A similar feature was found in 25 1 
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FIG. 3: (color online). HHG spectra calculated with the hard- 
core potential The black line represents the spectrum for 
muonic hydrogen at the laser parameters 7 H = 1.05 x 10 23 
W/cm 2 and fkj H — 59 eV. The grey (red) line represents the 
spectrum for muonic deuterium at the correspondingly scaled 
laser parameters I D = 1.30 x 10 23 W/cm 2 and hup = 62 eV. 
The inset shows a blow-up of the cutoff region on a linear 
scale, with the black line dashed for better visibility. 



and attributed to the behavior of non-tunneling harmon- 
ics in very steep potentials. More important with regard 
to the present study is, however, that the harmonic sig- 
nal from muonic hydrogen is larger (by about 50 % in 
the cutoff region) than that from deuterium. The reason 
is that a smaller nuclear radius increases the steepness 
of the potential near the origin, leading to more violent 
acceleration and thus to enhanced harmonic emission. 
With regard to ordinary atoms it has been noticed be- 
fore that the analytic behavior of the binding potential 
near the origin is of great importance for the HHG pro- 
cess [23j]; according to Ehrenfest's theorem, the atomic 
dipole is accelerated by the potential gradient which is 
largest in this region. In the present situation we there- 
fore find that the plateau height of HHG spectra from 
muonic atoms is sensitive to the finite nuclear size. 

In conclusion, we have found characteristic nuclear 
signatures in the harmonic response of strongly laser- 
driven muonic hydrogen and deuterium atoms. For deu- 
terium, due to the larger nuclear mass and size, both 
the harmonic cutoff and the plateau height are signif- 
icantly reduced. It was also shown that high-energy 
photons in the MeV regime can be generated which are 
capable of nuclear excitation. Our results could ex- 
perimentally be tested by virtue of near-future technol- 
ogy, i.e. by employing intense high-frequency radiation 
from free-electron lasers or plasma surface harmonics 
or, alternatively, by combining radiation from existing 
high-power near-optical laser systems with a relativis- 
tic atomic beam. The present study demonstrates that 
muonic atoms in strong laser fields can be utilized to dy- 
namically gain structure information on nuclear ground 
states. It also offers the prospect of performing pump- 



probe experiments on excited nuclear levels since the pe- 
riodically driven muon can excite the nucleus during one 
of the encounters and afterwards probe the excited state 
and its deexcitation mechanism. 
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